On non-Markovian time evolution in open 
quantum systems 

Andrzej Kossakowski 

Institute of Physics 
Nicholaus Copernicus University 
87-100 Torun, Poland 
kossak@fizyka.umk.pl 

and 

Rolando Rebolledo * 

Laboratorio de Andlisis Estocdtico 
Facultad de Matemdticas 
Pontificia Universidad Catolica de Chile 
Casilla 306, Santiago 22, Chile 
rreboUe&puc.cl 

Abstract 

Non-Markovian reduced dynamics of an open system is investigated. 
In the case the initial state of the reservoir is the vacuum state, an ap- 
proximation is introduced which makes possible to construct a reduced 
dynamics which is completely positive. 

1 Introduction 

An open system is one coupled to an external environment [U [2] ■ Such systems 
are of fundamental interest since the notion of an isolated system is almost 
always an idealization and approximation. The interaction between the system 
and its environment leads to phenomena of decoherence and dissipation, and 
for this reason recently received intense consideration in quantum information, 
where decoherence is viewed as a fundamental obstacle to the construction of 
quantum information processors [3]. 

In principle, the von Neumann equation for the total density matrix of the 
system and the reservoir provides complete theoretical predictions for all the 
observables. However, this equation is impossible to solve in practice because 
it takes into account all degrees of freedom of the reservoir. Efforts have been 
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focused on developing direct methods for the reduced density matrix of the 
subsystem. 

Two well known exact theories of subsystem dynamics are the Feynman- 
Vernon influence functional theory [H [51 [5] and the Nakajima-Zwanzig master 
equation [71 [5]. 

The Feynman- Vernon theory expresses the time evolving reduced density 
matrix of the subsystem as a path integral over subsystem trajectories weighted 
by an influence functional which incorporates the effects of the reservoir. In 
order to calculate the influence functional a path integral over all the reservoir 
degrees of freedom have to be performed. 

The Nakajima-Zwanzig approach employs projection operator techniques to 
derive an exact equation for the reduced density matrix from the von Neumann 
equation for the total density matrix. The resulting master equation -an integro- 
differential equation- is mostly of formal interest since such an exact equation 
can almost never be solved analytically or even numerically. In contrast, when 
one makes the Markovian approximation, i.e. when one neglects all reservoir 
memory effects, the resulting master equation [Sl[in] is formally solvable. More- 
over, the required property of complete positivity [11] is maintained. A coverted 
goal of the theory of open quantum systems is a non-Markovian description of 
time evolution which could at the same time include reservoir memory effects, 
remaining analytically tractable and retaining complete positivity. 

A variety of non-Markovian master equations have been proposed (cf. [1, 12 ~ 
32]). However, the complete positivity of the resulting time evolution is still an 
important problem to be investigated. 

On the other hand, in atomic, molecular and nuclear physics one deals with 
perturbation of discrete energy levels embedded in continuous spectra. Gamov 
[55] conjectured the existence of eigenvectors of the Hamiltonian corresponding 
to complex eigenvalues. However, it is impossible within the Hilbert space 
formulation of Quantum Mechanics because the Hamiltonian — being a self- 
adjoint operator — can have only real spectrum. In the case of simple scattering 
resonances Bohm and Gadella constructed the corresponding Gamov vectors 
in terms of a suitable extension of quantum theory on the basis of Gel'fand 
triples (see [35] [35]). This result was a first step towards a rigorous treatment 
of irreversibility [37-45]. Physical and mathematical aspects of Gamov states 
are presented in detail in [46| . 

In the present paper an attemp is made to apply these concepts to reduced 
dynamics of an open system. It is based on the observation that the spectrum 
of the free Hamiltonian of the system is embedded in the continuous spectrum 
of the reservoir and as a result of the interaction, the spectrum of the system 
becomes unstable, which leads to an irreversible evolution of the system. 

2 Reduced Dynamics in the Heisenberg Picture 

Let us consider a finite quantum system S with underlying complex separable 
Hilbert space f)"^. The reservoir R will be taken as an infinite quantum system 
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with the Hilbert space f)^ determined by the GNS representation tt^ induced by 
the reference state to, which is assumed to be invariant under the free evolution 
of R. 

The composed system S + R is considered to be isolated, i.e., its time evo- 
lution is determined by a bounded self-adjoint Hamiltonian Hx defined on the 
space f)'^ (g) f)^: 

Hx = (g) 1^ + l'^ (g, H"- + XV = Ho + XV. (1) 

The GNS representation associates to the state uj a vector in t)^ which we denote 
by the symbol \uj). Let P^^ = \iLj){i^\ be the projection on the state ui G 1)^. In 
terms of P^i, we define the following two projectors on the total space f)'^ (g f)^: 

Po = I'^f^P^ (2) 
Pi = 1-Po- (3) 

Notice that the following relation is satisfied: 

P^Ho = HoPo., (4) 

for a = 0, 1. 

Moreover, it is assumed that 

H'^Ilo) = 0. (5) 

The reduced dynamics % : £(f)s) — > ■C(f)s) and %t '■ 3i(t)s) 3i(f)s) in the 
Heisenberg and the Schrodinger picture, respectively, is defined by the relations: 

{ip,Ttia)^) = (Ut^(g,u,ia(g,l^)Ut^®Lo) = tr (ar,t(|7/;)((^|)) , (6) 

for all ip^tp & i)s and any a e £(f)s)- where 

Ut = eM^^tHx), (7) 
\^t) = Ut{\ip)®\u;)), (8) 
IV't) = (9) 

equation ^ can be rewritten in the form 

= {Po^t,{a®'i-'^)PoA) + {PiVt,{a'E)l''')Pii^t). (10) 

It follows from the above relation that the reduced dynamics is completely 
determined by the solution of the Schrodinger equation 



dt 

with the initial condition 



- -^Hx\^.), (11) 



lim|(^t) = \^)<»\uj). (12) 
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Using the projectors Pq, Pi and ^ one finds that and (|12p are equivalent 
to the following system of differential equations 

dPolv't) ^ tP„H^p„Po\^t)-iXPoVPiPi\^t), (13) 



dt 
dPiWt) 
dt 

with the initial conditions 



-iXPiVPoPoWt) - PiHxPiPiln), (14) 



frnPol^t) = (15) 
limPi|(^t} = 0. (16) 

The "variation of constants" method applied to (|14p with (1161) gives -Pi|iy9t) 
in terms of Po\(pt)- 

t 

Pi\iPt) = -*A J (e-'^'-^^''^"-P^PiVPoPo\^s}) ds, (17) 



which inserted in (|13p yields 



(18) 

Po|(^f} is first obtained solving (fTS]) under the initial condition (fT5| . After 
that, Pi|(pt) follows from (fT7|) and consequently, the reduced dynamics (fTO|) is 
determined. This dynamics is completely positive by definition. 
Notice that Po\ipt) is determined by 

PoUtPo = Poe-"«^Po, (19) 

and ([T8|) can be rewritten in the form 

t 

= -^PoH,PoPoU,Po-X' J (PoFPie-'(*-^)^^^^^^PiFPoPoC/.Po) ds . 



(20) 

Let denote PoU{p)Po the Laplace transform of PoUfPo, that is, 

PoU{p)Po = Poip + iHx)-'Po. (21) 

This means that PoU{p)Pq is the reduced resolvent. On the other hand, taking 
the Laplace transform of ([20| and using ((2T|) it follows that 

^o(P + ii^A)"'^0 = [p + iPoi/APo + A2PoyPl(p + lPll^Pl)-^PiyPo]"'Po, 

(22) 

which can also be derived via the resolvent equation. 

It is clear that the properties of Po\(pt) and Pil^pt) are determined by ana- 
lytical properties of the reduced resolvent. 
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3 The Friedrichs Approximation 



It is easy to understand that the formahsm developed so far can not be apphed in 
practice due to the comphcated nature of the unitary operator exp(—iPiHxPi) 
which appears in pT|) . and (|20p . In fact, the memory terms in p?)) . contain 
an infinite number of muhi-time correlation functions. 

It is worth noticing that the expression becomes simpler as soon as one 
assumes the additional condition 

PiVPi = 0. (23) 

Condition (|23p will be called the Friedrichs condition since it is satisfied in 
the Friedrichs model |47J as well as in its iV-level versions (see [48-52]). The 
simplification is due to the fact that in this case, 

PiHxPi - PiHoPi - PiHo. (24) 

The meaning of condition (|23p is that only one subspace of ()^ plays an 
essential role. 

In general the interaction Hamiltonian V does not satisfy (|23p. but V can 
always be decomposed as follows 

V = V ~ PiVPi+ PiVPi = V + PiVPi, (25) 

where V ^V- PiVPi ^ PqV + VPq - PaVPo- This operator has the following 
properties: 

PiVPi = 0, (26) 
Pi^Po = PiFPo, (27) 
PoVPi = PoVPi. (28) 

So that the Friedrichs approximation consists of replacing the interaction Hamil- 
tonian V hy V. 

In the case f)"^ — C", the interaction Hamiltonian may be written in the 
form: 

n 

hj = l 

where |ei), . . . , |e„) is an orthonormal basis in C" and the operators Fij satisfy 
Fij* = Fj,i. The operator V is then expressed by 

n 

that means that the only transitions in the reservoir are 
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To illustrate how the Friedrichs approximation works, let us consider an n- 
level system interacting with an electromagnetic field in the vacuum state. The 
corresponding Hilbert space for S is f)^ = C" and the Hamiltonian can be 
written in the form 

n 

= Y,^i\ei){ei\. (31) 

i=l 

For simplicity, it is assumed that the spectrum of is non-degenerate. 

The reservoir is chosen to consists of the quantized electromagnetic field. The 
modes of the electromagnetic field are indexed by A: = (k. A), where k e R^, 
\ € Z2 = {1, 2}. The Hilbert space of one photon states is f)i = L^(M^ ® Z2), 
and for any element / € [)i we write 



/ f{k)dk = V / /(k,A)d3k. 



At zero temperature, the Hilbert space of pure states of the reservoir is the 
Fock space ^ = ^k, where f)o = C and i)k is the symmetrized fc-fold tensor 

product of In 5" the creation and annihilation operators are introduced in 
the standard manner. In particular, 

a{f) = J a{k)f{k)dk, (32) 

for any f £ i)i and the vacuum state is determined by the condition 

a{f)\Lo) = 0, (33) 

for all / G f)i. 

The free evolution of the reservoir is defined by the formal Hamiltonian 

= j \k\a*{k)a{k)dk. (34) 
The simplest interaction Hamiltonian V can be chosen in the form 

n 

where fij S 

Making the Friedrichs approximation, i.e. V t-^ V — PqV + VPq — PqVPo 
with 

Po = 1^ (^\uj){uj\, 

one obtains 

n 

^ = Jl {\e,){ej\^a*{fi,j)\u){w\ + \ej){ei\ {u\a{fij)) . (36) 
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Now, applying (HOI) or one finds 



(37) 



where 

(e.,G(p)e.) = 5,AP + ^^k) + X'f: E / ■^'"■^■^^'.^^^-'-^f^' d^k . (38) 

The properties of G(p) depend on the choice of the functions /fe^j»(k, A). The 
above model is an analogue of the one due to Friedrichs. 

It should be pointed out that the solvability of the above model is related to 
the invariance of the subspace C" O w © C" ® [}i under the action of Hq + V. In 
spin-boson model (c.f. [2]) as well as models considered in [52-57], the existence 
of such an invariant subspace is due to an additional constant of motion. 



4 The Friedrich Approximation and non-Markovian 
Master Equation 



From (|10p it follows that the reduced dynamics in the Heisenberg picture can 
also be written in the form 

Ttia)®P^ = Poe''"ia^l'^')e-'*"Pa = Poait)Po, (39) 

where a(t) is the solution of the Heisenberg equation 

da{t) 
dt 

with the initial condition 



^ i[H,a{t)], (40) 



lima(t) = a®l^. (41) 

It is interesting to find out the master equation for 7^ (a) (8) P^j under the hy- 
pothesis that the interaction Hamiltonian satisfies the Friedrichs condition of the 
previous section, that is V has the form V = PqVPq + PqVPi +PiVP[j. For any 
X G £(fi), we denote Xc^p = P^XP^ for a, /3 = 0, 1. So that V = Voo+Vbi-hVlo- 
Thus, (PP]) becomes equivalent to the following system of equations: 



dt 
daoi{t) 

dt 
dawit) 

dt 
daiijt) 

dt 



= i[Hoo,aoQ{t)] + iiVoiaioit) - aoi(t)Fio) , (42) 

= i{Voiaii{t) - aoo{t)Voi) , (43) 

= i{Vioaoo{t) ~ an{t)Vio) , (44) 

= i[H°^,aii{t)]+i{Vwaoiit)~aw{t)Voi), (45) 
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with the initial conditions 



limaoo(t) 
lim aoi(t) 



hmaio(t) = hmaii(t) 



Ehminating aio(t) and aoi(t) in the previous system yields 



<Paii{t) 



with the initial copnditions 



00, 


dam{t)' 


dt \ 


rO 


daii{t)- 


111 


dt \ 



lim aoo(t) 

lim^ 

t— dt 

lim aii(t) 



t^o dt 

Finally, using the following matrix notations: 
\ aii{t) J ' 



+ 2V^)iaii(t)Fio - {aooit), VqiVw} 
+ 2T/ioaooWVbi - {aii(i), VwVoi} 

= a(g) , 

= 0, 
= 0. 



(46) 
(47) 



(48) 
(49) 

(50) 

(51) 
(52) 

(53) 



V = 





the above system of equations can be expressed in the form 
d^Ait) dA{t) 



Vol 
Vio 



= V* 



A^[V,[V,A(t)]]. 



dt^ dt 
Equation ([54|) is local in time but rather untractable. 



(54) 
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